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We quantify the capability of creating entanglement for a 
general physical interaction acting on two qubits. We give a 
procedure for optimizing the generation of entanglement. We 
also show that a Hamiltonian can create more entanglement 
if one uses auxiliary systems. 
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In the last 2-3 years there has been very considerable 
increase in experimental activity aiming to create entan- 
gled quantum states. One reason is the potential appli- 
cations of entanglement to quantum information process- 
ing. Creating entanglement has been possible in quan- 
tum optics for more than a decade, however now many 
new communities, working in a variety of experimental 
areas (for example NMR, condensed matter physics) are 
also joining the field jij. In general, entanglement be- 
tween two systems can be generated if they interact in 
a controlled way. However, in most experiments these 
interactions are weak which make the production of en- 
tanglement a very difficult task. Thus, it would be very 
convenient to have a theory which would provide us with 
the best way of exploiting interactions to produce entan- 
glement. 

In this Letter we analyze the entanglement capabilities 
of Hamiltonians. In particular, we would like to answer 
questions like: given an interaction (Hamiltonian), what 
is the most efficient way of entangling particles? Can 
we make the process more efficient by supplementing the 
action of the Hamiltonian with some local unitary opera- 
tions? Can wc increase the entanglement more efficiently 
by using some ancillas? 

So far, much of the theoretical effort in quantum infor- 
mation theory has been devoted to the characterization 
and quantification of the entanglement of a given state. 
Very recently, it has been realised that there is a parallel 
notion of the entanglement in the dynamics of a system 
In the authors consider the situation that one 
has a given unitary transformation and ask, for exam- 
ple, how much state entanglement is needed to produce 
it. Here we focus on a different issue: given an inter- 
action (i.e. a Hamiltonian) how can we make the most 
effective use of it ||. What we propose here is to define 
and determine the entanglement capabilities of physical 
processes, in particular, of unitary evolutions This 
is a very relevant problem not only from the theoretical 



point of view, but also from the experimental one. Of 
course, this problem is even more difficult than the one 
of quantifying the entanglement of states. In any case, in 
the present work we give the first steps in this direction 
by considering the case in which the physical process is 
acting on two qubits. 

From our results it turns out that: (i) It is more effi- 
cient to produce entanglement if one initially has already 
some; (ii) The best initial entanglement is universal, i.e. 
independent of the physical process; (iii) One can im- 
prove the performance of a physical process by comple- 
menting it with fast local operations; (iv) One can also 
improve it (in certain cases) by using auxiliary systems; 
(v) All entangling Hamiltonians can simulate each other 
and are thus qualitatively equivalent; we also provide an 
upper bound on the time required for one Hamiltonian 
to simulate another. 

We consider two qubits interacting via a non-local 
Hamiltonian H . We want to determine the most efficient 
way in which we can use such an interaction to produce 
entanglement. We will characterize the entanglement of 
a state of the qubits at a given time t, \9(t)), by some 
entanglement measure E. In order to quantify the en- 
tanglement production, we define the entanglement rate 
r at a particular time t of the interaction as follows: 



r(t) 



dE{t) 
alt 
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This quantity depends on | not only through its 
entanglement E. The goal is then to find the conditions 
which must be satisfied in order to obtain a maximal 
entanglement rate. In particular, we will be interested in 
determining the following: 

(i) For any initial entanglement E of the two-qubit sys- 
tem, what is the state | , F), say \^e), for which the 
interaction produces the maximal rate Te- 



(ii) The maximal achievable entanglement rate r max [pi , 
r max = max^T^ and the state |\l/ m ax) for which 
F — F 

i x max ■ 



These quantities are interesting because the knowledge 
of the state \^e) will allow us to find out the most effi- 
cient way of entangling the qubits. The idea is to sup- 
plement the interaction Hamiltonian H with appropriate 
local unitary operations in such a way that the state of 
the qubits at any time t is precisely \*&E(t)}, f° r which the 
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increase of entanglement is optimal. In order to show how 
this can be achieved, let us consider that the evolution 
given by H proceeds in very small time steps St. Let 
us also assume that the qubits are initially disentangled. 
Using local operations, we can always prepare the state 
|\&o) — that is, the product state which most efficiently 
becomes entangled under the action of H. After a time 
step St, the state will change and its entanglement will 
increase to SE. Then, we use (fast) local unitary oper- 
ations to transform the new state of the qubits into the 
state I^e) for which T is optimal. Note that this is al- 
ways possible, since for qubits all states with the same 
value of E, say SE, are connected by local unitary trans- 
formations. By proceeding in the same way after every 
time step, and taking the continuous time limit <5t — > 
we obtain that the state of the qubits at time t is always 
the optimal one, \$f E(t))- Obviously, in an experimen- 
tal realization, this procedure requires that we can apply 
the appropriate local transformations in times which are 
short compared to the typical time scale tjj associated 
to H, r H = (e max — emin) -1 , where 

e max and e m i n are the 
maximum and minimum eigenvalues of H, and we have 
set h = 1. 

Knowledge of Te also permits us to determine the 
maximum amount of entanglement E max produced as a 
function of time. We just have to express as an ex- 
plicit function of E, substitute it in (|l|) and solve that dif- 
ferential equation to determine E max (t). Note that the 
optimal procedure described above will precisely reach 
the entanglement E max (t). 

The state l^max) is important since it gives rise to the 
maximal increase of entanglement, and therefore corre- 
sponds to the best operational point. After reaching the 
state l^max) with the procedure described in the previ- 
ous paragraph, the entanglement would be produced in a 
very efficient way, if one could transfer the entanglement 
that is gained after each time step St to other qubits (us- 
ing entanglement dilution || or some other means). In 
particular, it would increase proportionally to the time, 
r max being the proportionality constant. 

In the following, we will show how to determine I^b), 
Tb, l^max), and r max for an arbitrary Hamiltonian H. 
To this end, it is convenient to use the Schmidt decom- 
position of the state of the qubits to write 



P\<p, X )+e* a VT^P\<p ± ,x ± ), 



(2) 



where for the sake of short-hand notation we have omit- 
ted the time dependence of all these quantities. Here, 
(ipl^) = (x\x X ) = and P < 1/2. Note that E must 
only depend on the Schmidt coefficient P, given the fact 
that it must be invariant under local unitary operations. 
For example, if we choose as entanglement measure the 
entropy of entanglement || — the entropy of the reduced 
density operator of one of the qubits — , we will have 



E(P) = -Plog 2 (P) - (1 - P) log 2 (l - P). 



(3) 



Note that the entropy of entanglement quantifies the 
amount of EPR entanglement contained asymptotically 



in a pure state \^f). That is E(P) gives the ratio of max- 
imally entangled EPR states |*") = l/v^flOl) - |10)) 
which can be distilled from [are needed to create] \^>) 
respectively. Thus, we can write 



W dP dt 



(4) 



In (Q), given a particular entanglement measure E(P), 
we just have to determine dP/dt. In order to do 
that, we need to find the (infinitesimal) time evolu- 
tion of the Schmidt coefficients of the state of the 
qubits. After a time St we will have \^f(t + St)) = 
exp(-iH5t)\^(t)) ~ (1 - iHSt)\^{t)). The corre- 
sponding reduced density operator pA{t + St), where 
Pa,b — ^b,a{\^)(^\), can then be written as pA{t + 
St) = PA (t) - iStTr B {[H, \$>(t))(^(t)\}}. The eigenvalues 
(Schmidt coefficients) of this operator can be easily de- 
termined starting from pa\<p) = PW) and using standard 
perturbation theory. We find 



dP 
~dt 



= 2^/P(l - P) x Im[e ia {<p, xWl? , X^} 



(5) 



where we have omitted the time-dependence. Upon sub- 
stitution in (Q) we obtain the entanglement rate. Since 
we are interested in maximizing T, it is clear that we can 
always choose a such that 



T = f(P)\h(H,cp, X )\. 



where 



f(P) = 2^P(1-P)E'(P), 
h{H, V , X ) = (<p,x\H\ V X ,X X )- 



(6) 



(7a) 
(7b) 



By analyzing Eq. (|g) we can extract some interesting 
conclusions, even before determining the maximum value 
of r explicitly. Given the fact that / and h depend on 
different parameters, in order to determine the quantities 
mentioned in (i-ii) we can maximize the functions / and 
\h\ independently. First, if we want to determine the 
quantities mentioned in (i), we have to fix the value of E. 
In that case, P is also fixed and therefore the maximum 
of the entanglement rate will correspond to a state of the 
form (^|) with some fixed | tp) , \x), and a (which maximize 
\h\). That is, for any value E of the entanglement, the 
states \ip) and |x) for which the maximal entanglement 
rate is obtained do not depend on E, but only on the 
form of the Hamiltonian H. Let us denote by /i max the 
maximum value of \h\; that is, 



IMI, llxll=i 



Then, we can easily determine how the entanglement 
would evolve with time if we always drive the qubits with 
local operations so that at each time their state corre- 
sponds to the optimal one. We can simply solve the dif- 
ferential equation (S), obtaining P{t) = sin 2 [/i max t + O ], 
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with P(0) = sin 2 



Using the explicit dependence of 



E on P, we can directly then calculate E(t). The evolu- 
tion of the entanglement is fully characterized by /i max , 
which is a quantity that only depends on the interaction 
Hamiltonian. That is, for a given H, /i max measures the 
capability of creating entanglement. In the following we 
give a simple way of determining ft. max , which allows us 
to classify the entanglement capability of any Hamilto- 
nian. On the other hand, once the entanglement mea- 
sure is specified, we can calculate the value Po of P for 
which we obtain the maximal rate by simply considering 
the function f(P). For example, choosing the expres- 
sion (|^) for the entanglement, we find that Po solves the 
equation In 1 p^ >0 = 1 _2 Pg > i- e - Po — 0.0832 which gives 
P(Po) ~ 0.413. This shows that, in order to increase the 
entanglement of a two-qubit system in an optimal way, 
it is better to start with some initially entangled state 
rather than a product state 0. Note that the optimal 
initial entanglement E(Pq) is independent of H. 

In the following we will show how to determine the 
entanglement capability /i max of a general Hamiltonian 
H acting on the qubits. First, we will show how, if wc 
allow to supplement the evolution of H by local unitary 
operations, we can express H in a standard form that 
only depends on three parameters. Then we will derive 
an expression for /i max in terms of those parameters. 

Except for a trivial constant, we can always express a 
general Hamiltonian as 

3 3 

H = &iot ® 1 b + ^ Po 1 ^ ® of + 

i=l j=l 



J2 7i,j°f®° 



(9) 



Here, Oi are the Pauli operators, and a, (3, and 7 are two 
real vectors and a real matrix, respectively. We now show 
that by supplementing the evolution operator with local 
unitary operations we can obtain an effective Hamilto- 
nian which has one of the two standard forms 

H ± = fnaf <g> erf ± ® erf + fi 3 v£ ® erf , (10) 



3 > are the (sorted) singular values 
1. We first note that the terms cor- 



where /ii > /12 > 
of the matrix 7 
responding to a, f3 in @ give no contribution to /i max 
(||) and can therefore be neglected. Second, we apply 
the local operations U (V) and W (V^) to the first 
(second) qubit at the beginning and end of the evolu- 
tion process, respectively. We select them such that 

u^fu = Zl =1 o£X> v ^f v = Ef=i°jK< where 

qA.b are or thogonal matrices of determinant one, each 
being plus or minus the orthogonal matrices in a singular 
value decomposition of 7. Thus the total (non-local) ef- 
fect of the evolution for a time t is equivalent to the one 
obtained with the Hamiltonian H + (H~) for the same 
time if det7 > (det7 < 0) . Without loss of generality, 
we may take H of the form H+ (jic]) in what follows §. 



Now let us determine /i max in terms of /ii,2,3- We can 

write h(H,ip,x) = Efe=i l*k{f |of Iv^ixWk \x X )- Using 
the Cauchy-Schwarz inequality, it can be checked that 
the maximum of (the absolute value of) this function is 
reached for \x) = \(f )■ In this case, using the fact that 
|V3)((^| + |^- L )(<^- L | = 1 we obtain h(H, ip, ip) = Y^l^i^- 

Efc=i A t fc(v 5 l cr fc|'/ 9 ) 2 - Taking into account that /ii > ^2 > 
/X3, we see that the maximum value occurs when \tp) = |0) 
or = |1), i.e. an eigenstate of 03. For that choice we 
obtain 



/imax — Mi + f-2 ■ 



(11) 



Summarizing, once we have transformed the Hamiltonian 
H to the standard form (|l(]) we obtain that for a given 
value of E (and therefore of P), 



/ P|0,l)+i v / l _ P|1,0), 

/(P)<Wx 



(12a) 
(12b) 



where /i max = fix + /i2- The maximum rate r max is ob- 
tained for P = Pq, where Po is the value that maximizes 
/(P). Thus, |* m ax) and r max are given by (|l|) with 
P = Pq. For example, for the entanglement measure (|^), 
P ~ 0.0832 which leads to /(P ) ~ 1.9123. 

So far, we have calculated the most efficient way of 
entangling two qubits if we can use local unitary oper- 
ations acting on each of the qubits. We have not al- 
lowed, however, local operations which entangle each of 
the qubits with local ancillas. We will now show that 
this possibility permits us to increase the maximum en- 
tanglement rate r max for certain kind of Hamiltonians. 
We will first generalize the formulas derived above to 
the case of multilevel systems, given that the system 
qubit-plus-ancilla is of this sort. We consider a state 
with Schmidt decomposition l^) = J2n=i V^nl^n, Xn)- 
As before, any entanglement measure E will only de- 
pend on the Schmidt coefficients A„ > 0. In particular, 
in the following we will use the entropy of entanglement, 
E(X) 

entanglement rate, we have 



J2n=i log 2 (A„). Using the definition (0) of 



N 



dE d\ ri 



1 



^-f dX n dt 

n— 1 



N 

N ^ 

n.m— 1 



dE 
dX n 



dE 



dX„ 



dX n 
dt ] 



(13) 



where we have used the fact that the sum of all the 
Schmidt coefficients is constant. Using perturbation the- 
ory as before, we find ^f- = 2V%Im[(<£„, x«|P|*>] = 

2 Em=l V*>nX m Im[(ip n , Xn\H\ip m , Xm)]- 

Rather than proceeding in complete generality we now 
consider an example which demonstrates that adding an- 
cillas may allow one to increase entanglement more effi- 
ciently than is possible without the use of ancillas. We 
will consider the case in which the ancillas are also qubits. 
We write P — X\ and concentrate on the case in which 
A2 = A3 = A4 = (1 — P)/3. In that case, Eq. ( |i~3| ) simpli- 
fies to 
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f = f(P)h(H,cp n , X n) 



(14) 



where now 



/(P) = 2y/P[l - P)/31og 2 [(l - P)/(3P)], (15s 



h(H,(f n ,Xn) = ^2lm[(ifii, Xl\H\ifi n , Xn)}- 



(15b) 



n=2 



We can always choose the phase of the states \ip n ) such 
that all the terms on the sum add with the same sign. Wc 
can therefore replace the imaginary parts of the terms 
in the above expression by their absolute values, and 
in (|J) we can replace /(P) by |/(P)|. Wc find that 
Po ~ 0.8515 (which corresponds to an entropy of entan- 
glement E(P ) ~ 0.8415) maximizes \f{P)\ ( p^aj ) and 
leads to |/(Po)| — 1.6853. Proceeding as before, we can 
easily maximize h. We obtain that the maximum value 
is hmax = Mi + M2 + 1*3, which occurs when \(p n ) = \x n ) 
are orthogonal maximally entangled states between the 
qubit and the ancilla. For example in the case that 
det7 > 0, the choice \ipi) = \(f> + ), \(f2) = l^s) = 

ji Up~), \ f4) = where {|0 ± ), |V' ± )} are Bell states 

jUjl , together with P = Pq = 0.8515 leads to a maximal 
[under the previous assumptions on the Aj's] entangle- 
ment rate T = r max . 

Let us compare the cases in which we use ancillas 
and the one in which we do not use them. On the one 
hand, we have that |/(Po)| < |/(Po)|. But on the other, 
/i max > /i m ax- Thus, if /i3 7^ it may be the case that the 
use of ancillas can help to increase the maximum rate of 
entanglement r max as well as the rate Te for a given en- 
tanglement E of state This is in fact the case, if we 
have, for example, ^1=^2 = (i-e. /i max = 3/2/i max ). 
In this case we obtain r max ~ 1.3220r max . In a simi- 
lar way, one can check for this specific Hamiltonian that 
Tb > r_E if the initial entanglement satisfies E > 0.08. 

Finally, it is easy to show that all entangling Hamilto- 
nians are qualitatively equivalent when assisted by local 
operations. In particular, given two Hamiltonians H and 
W with either /i max = ah' max or ft, max = ah' max , one can 
simulate the action of H' for any time t by applying H 
for at most 3a -1 £. This can be seen as follows: Applying 
H of the form (|l0|) for St/2 followed by a local unitary 
operation <j\ in A before and after another application of 
H for St/ 2 is equivalent to the application of the Hamil- 
tonian Hi = [L\0\ ® o\ for the time St, provided that St 
is infinitessimally small. Since Hi is locally equivalent to 
Hk = MiCfc^c/c, applying sequentially Hk for 5tfj,' k /fXi is 
equivalent to the application of H' for the time St. Using 
the restrictions on ^fc,/Lt' fc the claim readily follows. The 
question of efficient simulation of another Hamiltonian in 
the same time t will be addressed in future work. 

In summary, we have found the optimal way of using 
any non-local interaction to entangle a pair of qubits. 
The idea is to use local operators to drive the instanta- 
neous state to the one that maximizes the entanglement 



rate, at each moment of the evolution. We have found 
that the entanglement capacity of any given Hamilto- 
nian is determined by the sum of the two largest singular 
values of the matrix 7 defined in (^). Finally, we have 
shown that for certain Hamiltonians one can overcome 
this maximal entanglement rate by using ancillas pre- 
pared in maximally entangled states with the qubits. 
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Note added: The idea of universal simulation of Hamil- 
tonians as discussed in this manuscript originated in a 
larger collaboration, including several other people and 
is developed in detail in p| . After completion of this 
manuscript we have also learned that the notion of uni- 
versal simulation of Hamiltonians has been independetly 
addressed by Dodd et al. |l2| . 
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